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We investigate hadronic unquenching effects in light quarks and mesons. Within the non- 
perturbative continuum framework of Schwinger-Dyson and Bethe-Salpeter equations we quantify 
the strength of the back reaction of the pion onto the quark-gluon interaction. To this end we add 
a Yang-Mills part of the interaction such that unquenched lattice results for various current quark 
masses are reproduced. We find considerable effects in the quark mass function at low momenta as 
well as for the chiral condensate. The quark wave function is less affected. The Gell-Mann-Oakes- 
. Renner relation is valid to good accuracy up to pion masses of 400-500 MeV. As a byproduct of our 

investigation we verify the Coleman theorem, that chiral symmetry cannot be broken spontaneously 
f*"^ ' when QCD is reduced to 1+1 dimensions. 
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I. INTRODUCTION 



Dynamical chiral symmetry breaking is one of the most important properties of low- 
Q-f energy QCD. The breaking pattern has profound impact for phenomenological quantities, 
Q_i' as e.g. the appearance of the pseudoscalar Goldstone bosons in the chiral limit of QCD 
and the non-degeneracy of chiral partners. Chiral perturbation theory [l], [2| describes 
these effects very efficiently on the level of hadrons but has nothing to say about the 
\ underlying structure of the full theory. In this work we investigate dynamical chiral 
>- ■ symmetry breaking on the level of (non-perturbative) correlation functions between quarks 
and gluons. 

We are particularly interested in the interplay between the fundamental quark and 
gluon degrees of freedom and the resulting bound states. In full QCD there are hadronic 
lO ■ contributions to the fully dressed quark-gluon interaction. These effects are generated by 
the inclusion of dynamical sea quarks in the quark-gluon interaction and are therefore 
only present in unquenched QCD. In this work we concentrate on the back reaction of the 
pions on the quark propagator and investigate the impact of the unquenched interaction 
on the quark-pion system. 

Pion effects on the quark propagation are important for several reasons. They account 
for (at least part of the) pion cloud effects in baryons and mesons. Thus they need to 
be incorporated in bound state calculations aiming at an adequate description of phe- 
nomenological properties of these objects. Furthermore they allow for the possibility of 
hadronic intermediate states in bound state calculations and therefore generate the finite 
width of meson spectral functions. On a more fundamental level Gribov argued that pion 
effects in the quark propagator may be responsible for quark confinement [3j. In addition, 
the inclusion of these effects in a continuum framework allows for an extrapolation of 
unquenched lattice results for large quark masses towards the physical up/down quark 
values. Finally, it will become apparent that Goldstone effects are essential in smaller 
dimensions. In two dimensions they prohibit the dynamical breaking of a continuous 
symmetry Q]. 

We study the pion back reaction on the quarks in the Green's function approach to 
Landau gauge QCD using Schwinger-Dyson and Bethe-Salpeter equations (SDE/BSE) jl, 
@, 0, II] • In the past these effects have been investigated to some extent within the N JL- 



2 



model [§, QiJ U|. Here we are extending these model studies to full QCD. Unquenching 
effects in the quark-antiquark system have also been investigated by Watson and Cassing 
within the SDE/BSE approach in [12j]. There a coupled system of equations for the Bethe- 
Salpeter kernel and the connected quark four-point function has been solved in a certain 
approximation, which neglects the back reaction of the hadronic resonances on the quark 
propagator. Our study here is therefore in some sense complementary with their work. 
It is also complementary to the investigations reported in jl3l . Hj], where unquenching 
effects in the gluon polarization have been considered. 

In order to make our results as concise as possible we use a combination of two methods. 
We devise a truncation scheme for the combined Schwinger-Dyson equations for the gluon, 
ghost and quark propagator as well as the quark-gluon vertex, which is combined with the 
corresponding Bethe-Salpeter equation for light mesons. The details of the quark-gluon 
interaction are fixed such that the available unquenched lattice results for the quark 
propagator [l5| are reproduced. This combination of methods allows us to study the 
relative impact of the pion contributions compared to the pure Yang-Mills part of the 
quark-gluon interaction. 

The paper is organized as follows. In subsection III Al we discuss the hadronic contribu- 
tions in the quark-gluon vertex and detail our approximation scheme of this vertex and 
the quark Schwinger-Dyson equation (SDE). In subsection III Bl we explain the extrac- 
tion of pion contributions in the quark-antiquark scattering kernel of the corresponding 
Bethe-Salpeter equation (BSE). The resulting interaction is supplemented by the familiar 
one-gluon interaction from the rainbow-ladder approximation. Both types of interac- 
tion together represent an approximation to the full unquenched quark-gluon interac- 
tion, which retains important properties of the mechanism of dynamical chiral symmetry 
breaking. A different perspective to this approximation in the framework of 2PI effective 
actions is discussed in subsection III CI In subsection III Dl we explain our procedure to fit 
the parameters of the one-gluon exchange part of the interaction such that results from 
unquenched lattice simulations at large quark masses are recovered. Some technical de- 
tails of our numerical method are presented in IIII Al In subsection IIIIBI we discuss our 
numerical results for the quark propagator and the pion at physical quark masses. We 
then investigate the chiral limit in subsection IIII CI We show that the Gell-Mann-Oakes- 
Renner relation is satisfied in our approach and give a value for the chiral condensate. In 
subsection IIII Dl a preliminary analysis for the analytical structure of the resulting quark- 
propagator is presented. Our results are summarized and discussed further in section IIVI 
Throughout the paper we work in Euclidean space and in the isospin symmetric limit of 
equal up and down quark masses. 



II. THE QUARK-GLUON INTERACTION 



As emphasized above there is a distinct imprint of dynamical chiral symmetry break- 
ing on the low-energy meson spectrum: the lowest-lying pseudoscalar meson states are 
identified with the (pseudo-)Goldstone bosons of this mechanism. In the chiral limit, 
these states are exactly massless although they can also be described as bound states of 
constituent quarks with masses of the order of 300-400 MeV. This dichotonuous nature of 
the pseudoscalar states has been discussed from first principles in Q . It is the axialvector 
Ward-Takahashi identity (axWTI), relating the quark self energy and the quark-antiquark 
scattering kernel, that enforces a binding energy of the pseudoscalar meson system which 
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exactly cancels the quark and antiquark masses. It is therefore mandatory for every 
meaningful computational scheme of QCD to respect the axWTI. 

A widely used practical truncation of the SDE/BSE framework with this property is the 
rainbow-ladder approximation for the quark interaction. It has a history of remarkable 
successes (summarized e.g. in [7]). However, there are also shortcomings that limit the 
credibility of such an approximation. Consequently, several efforts have been made to 



extend this scheme, see e.g. [l_4j, lla, [13, llSl, ll9|, |20|]. In the following we will extend the 



rainbow-ladder scheme by taking into account additional hadronic contributions to the 
quark-gluon interaction. 



A. SDEs for the quark-gluon vertex and the quark propagator 

Let us first briefly sketch the so-called renormalizat ion-group- improved rainbow-ladder 
truncation. The starting point is the SDE of the quark propagator 

s-\ P ) = z 2 s 1 ( p ) + j:(p) j (i) 

where Sq 1 (p) = ip-j + m denotes the inverse bare quark-propagator, while S~ 1 (p) = 
ip-^A(p 2 ) + B(p 2 ) is the dressed propagator and Z 2 the renormalization factor of the 
quark field. The quark self energies A(p 2 ) and B(p 2 ) can be recombined into the quark 
mass M(p 2 ) = B(p 2 )/A(p 2 ) and the quark wave function Zf(p 2 ) = 1/A(p 2 ). The quark 
self energy is given by 

S(p) = g 2 C F Z 1F JlA-^Siq^^^D^k), (2) 

with k = p — q, the Casimir Cf = {N 2 — 1)/(2N C ) and the renormalization factor Z\p of 
the quark gluon vertex. The self energy depends on the fully dressed quark-gluon vertex 
r„(<7, k) and the gluon propagator 

kjX\ Z(k 
/ '" Ah 1 \ !,/ "' k 2 J k 



DM = - , (3) 



with the gluon dressing function Zik 2 ). Throughout this paper we will work in the Landau 
gauge, which guarantees the transversality of (J3j). The rainbow-ladder approximation 
amounts to the replacement 

^Z(k 2 )T„(q,k)^^T(k 2 ) lu , (4) 

where T(k 2 ) can be viewed as a combination of the gluon dressing function and a purely 
/c 2 -dependent dressing of the 7„-part of the quark-gluon vertex. 

Often the function T(k 2 ) is referred to as an effective running coupling. Indeed it can be 
shown, that the ultraviolet behavior of T(k 2 ) has to resemble that of the strong running 
coupling in order to reproduce the correct ultraviolet behavior of the quark propagator 



as known from resummed perturbation theory and the operator product expansion [22 
The rainbow approximation then gives reliable results for the quark propagator in the 
ultraviolet. In the infrared, however, one needs to go beyond simple approximations 
like (jl]). The infrared shape of the quark propagator and also its analytical structure 
do depend on the details of the quark-gluon interaction [l3|, 23[. In particular, tensor 



structures other than 7^ seem to be important there [24 
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FIG. 1: The full, untruncated Schwinger-Dyson equation for the quark-gluon vertex [2l( is shown diagrammatically 
in the first line. The second line describes the first terms of an expansion in terms of hadronic and non-hadronic 
contributions to the quark-antiquark scattering kernel. In both equations, all internal propagators are fully 
dressed. Internal dashed lines with arrows correspond to ghost propagators, curly lines to gluons and full lines 
to quark propagators. In the second equation, the dotted line describes mesons, the dashed line baryons and the 
double lines correspond to diquarks. 



The Schwinger-Dyson equation for the quark-gluon vertex is given diagrammatically in 
Fig. [TJ In the first line we show the full, untruncated equation. For very small momenta, 



a selfconsistent solution to this equation has been given in ref. [241 ] . Here we are primarily 



interested in the mid-momentum behavior of the vertex and in particular in hadronic 
contributions. To lowest order in a skeleton expansion such contributions can only oc- 
cur in the diagram with the bare quark-gluon vertex at the external gluon line. In the 
second line of Fig. [TJ we expand the quark-antiquark scattering kernel of this diagram in 
terms of resonance contributions to the kernel and one-particle irreducible Green's func- 
tions ('skeleton expansion'). Amongst other terms one obtains one-meson and one-gluon 
exchange between the quark and anti-quark lines. Also diquark exchange contributions 
arise which are not explicitly shown in Fig. [TJ The first baryon exchange diagram shows 
up as a 'two-loop' diagram, which also involves diquarks and the baryon-quark-diquark 
vertices. Note that double counting is trivially avoided in this combined expansion due 
to different quantum numbers in the exchange channel. 

The computation of the hadronic diagrams in Fig. [TJ is rather involved. The meson- 
exchange diagram requires in addition the solution of a coupled system of the Schwinger- 
Dyson equation for the quark propagator and a corresponding Bethe-Salpeter equation 
(BSE) for the meson-quark vertex. Even more complex is the baryon-exchange diagram, 
which involves the computation of the diquark-BSE and a Faddeev-type equation for the 
baryon bound state. In this work we therefore wish to concentrate solely on contributions 
from the one-pion exchange diagram in addition to the non-hadronic contributions. In fact 
this choice should be regarded as a first approximation to the fully unquenched system, 
since all other hadronic diagrams are at least suppressed by factors of Ag^/m^ with 
H G {K, p, N, ...}. The same is true for diquark exchange contributions. 

On the quark level, the one-pion exchange diagram in Fig. [TJ involves a closed quark 
loop and therefore can only appear in unquenched QCD. Our investigation therefore 
complements previous studies of unquenching effects in the SDE framework, where quark- 



FIG. 2: The approximated Schwinger-Dyson equation for the quark-gluon vertex. All internal propagators are 
fully dressed. 
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FIG. 3: The Schwinger-Dyson equation for the quark propagator with the quark-gluon vertex from Fig. [2] 



loops in the gluon polarization have been investigated [13|, [14| . 



The Yang-Mills part of the vertex is conveniently approximated by 

r™(Pi, = i»r YM (pl), (5) 

where we denote the quark momenta by p\ and P2 and the gluon momentum by p^. 
The explicit expression for Y YM (p1) is detailed in subsection III Dt The ansatz (jSJ) is an 
example for the rainbow- ladder approximation, cf. Eq. (pO). A similar form has been used 



in |25l. l26l|. where quenched lattice results for the quark propagator have been reproduced. 
Note, however, that (j5J) involves only the 7^-part of the tensor structure of the full vertex. 

In general one can decompose the quark-gluon vertex Y u (pi,P2,P3) into twelve different 
tensor structures, given by 
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^u(Pl,P2,P3) =^2Xi(pi,P2,P3)K(PuP2,P3) ■ (6) 



i=l 



The details of this basis are given in ref. 27|. Here we only note that L\ = 7^; the explicit 



forms of the other structures are not needed. In principle, all twelve tensor structures 
can be important in the intermediate momentum regime. This has been explored to 
some extent in lattice simulations j28|. We have verified by projection methods that the 
one-pion exchange diagram of Fig. [1] contributes to all these structures. 

These considerations lead to the following approximation scheme for the quark-gluon 
vertex: we subsume the Yang-Mills part of the quark-gluon interaction into a form as 
given in Eq. (jSj) and add the one-pion exchange diagram of Fig. [TJ The resulting vertex 
is given in Fig. [2J It contains pion contributions in all twelve dressing functions Aj and 
contributions from the Yang- Mills sector in A x . A conceptually different justification 
of this approximation scheme in terms of a l/A^-expansion is given in subsection III CI 
Certainly, an explicit diagrammatic calculation of the Yang-Mills part of the vertex would 
be preferable to an ansatz of the form of Eq. (JSJ). Such a calculation is under way and 
will be detailed elsewhere. 

Inserting the vertex of Fig.[2]into the Schwinger-Dyson equation of the quark-propagator 
we arrive at the diagrammatic equation given in Fig. [31 In order to be able to construct a 
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FIG. 4: The Bethe-Salpeter equation in rainbow- ladder approximation, i.e. effective one-gluon exchange. 



corresponding kernel for the Bethe-Salpeter equation of the pion we need to perform an 
additional approximation as detailed in the next subsection. 

B. The quark-ant iquark scattering kernel 



We start with the general expression for the homogeneous Bethe-Salpeter equation 
(BSE) for quark-antiquark bound states, which can be written as 



d 4 k 



(2 



7T 



K tu . rs (p, k; P) [S(k + )T(k; P)S{k„)] s 



(7) 



Here T is the Bethe-Salpeter vertex function of a quark-antiquark bound state and K 
is the Bethe-Salpeter kernel. The momenta k + = k + £P and k- = k + (£ — 1)P are 
such that the total momentum is given by P = k + — k_. The momentum partitioning 
parameter £ reflects the arbitrariness in the relative momenta of the quark-antiquark 
pair and can be set to £ = 1/2 without loss of generality. Latin indices (t,u,r,s) refer 
to color, flavor and Dirac structure. The BSE is a parametric eigenvalue equation with 
discrete solutions P 2 = —M 2 where M n is the mass of the resonance. The lowest mass 
solution corresponds to the physical ground state. Since P 2 is negative, the momenta k± 
are necessarily complex in Euclidean space and thus the quark propagator functions are 
evaluated with complex argument. We will come back to this issue below. 

Chiral symmetry constrains the Bethe-Salpeter kernel K tu . rs via the axialvector Ward- 
Takahashi identity (axWTI), 



[E(p + ) T5 + 75 E(p_)] 



d 4 k 

W) 1 



K tu . sr (p, k; P) [j B S(k-) + S(k + )j 5 



rs ' 



(8) 



which relates the kernel to the quark self energy E(p) = S^ 1 (p) — Z 2 Sq X (p) , cf. Eq. (PQ). In 
the rainbow-ladder approximation this relation can be satisfied easily. For our Yang-Mills 
part of the interaction, Eq. ([5]), one then obtains the kernel 



K™(q,p;P) 



g 2 Z(k 2 )T YM (k 2 )Z 



IF 



k 2 



kp}z v 



"7m 



J ts 



2 lv 



(9) 



The resulting BSE includes an effective one-gluon exchange between the quark-antiquark 
pair as shown in Fig. HI 

The construction of a Bethe-Salpeter kernel corresponding to the pion-exchange part 
of our interaction is more complicated. We were not able to find an exact solution of the 
axWTI for the interaction given in Fig. [2j In principle one then has two choices: either one 
works with an approximate solution of the axWTI or one approximates the interaction 
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further such that an exact solution of the axWTI to the approximated interaction is 
possible. Since the first option destroys the Goldstone character of the pion, only the 
second choice is viable and shall be followed here. 

To this end we reconsider the quark-SDE given in Fig. [3j For diagrammatical reasons 
the corresponding Bethe-Salpeter kernel should look like a superposition of the effective 
one-gluon exchange and an interaction involving pions. Let us assume for the moment 
that only the Yang-Mills part of the interaction would be present in the BSE. Then the 
quark-SDE could be rewritten by inserting the BSE into the second diagram on the right 
hand side. The result is given diagrammatically in Fig. [5] and represents an interaction 
given by an effective one-gluon and a one-pion exchange diagram in the SDE. 

The only remaining question concerns the description of the one-pion exchange, which 
will be needed in the t-channel (see Fig. [3] and Fig. [5]). In the s-channel the pion contri- 
bution to the quark-antiquark scattering kernel is given by [B| 

M£?.{q,P\ P) = KU (?; -P) p^j^Khu (p; P) + O ((P 2 + M 2 ) ) , (10) 

— ' i T 

where is the Bethe-Salpeter vertex function, T l n (q; —P) = C Y\ (—q; —P) C and 
C = 7274. As already stated we neglect other contributions in the quark-antiquark scat- 
tering kernel as we expect the Goldstone modes to be dominant. The corresponding 
expression for the t-channel should be related by crossing symmetry (if we take incoming 
and outgoing particles on-shell). As the latter will however be broken when only consid- 
ering the leading term in Eq. ffTUl) . a transformation of the leading term in Eq. (fTUj) to the 
t-channel is not unique and depends on the choice for incoming and outgoing momentum. 
In order to find a suited kernel for the BSE which fulfills the axWTI in Eq. (jSj) later on, we 
take the arithmetic mean of the two possible contributions and obtain for the self-energy 

A cf , (, k^k u \ Z(k 2 )T YM (k 2 ) 



S 1 {p) = Z 2 S Q 1 {p)+g 2 C F Z lF J tJ^7/,<%)7^ f<W fc2 j fc2 



W) 



P + 1 \ Qf x r ( P + Q 

P-q] SiqjT^ — — ;q-p 



+ K [^—^;q-p] S(q)T 7T l^-^-;p- q 



D n (k 2 



with k = p — q and D n (k) = (k 2 + M 2 )~ x being the pion propagator. The factor 3 in front 
of the pion contribution stems from the flavor trace and represents contributions from 
7T + , 7r_ and 7r . These are treated on equal footing in the isospin-symmetric limit adopted 
here. The general structure of the Bethe-Salpeter vertex of the pion can be represented 
by 

Tl(p, P) = r l 75 {E n (p, P) - iPF^ip, P) - ip- 1P ■ PG n (p, P) - tfp.y] H n (p, P)} (12) 

with four independent dressing functions E n ,F n ,G n and H n . Those obey f(k;P) = 
f(—k; — P) = f(—k; P) for / e {E n , F w , G w , H T }. For the case H T = it is then possible 
to respect the axWTI in Eq. (JSJ) for P 2 — > in the chiral limit by choosing the kernel 

A-£i(9,p;-P) = ( P + q 2 +P ;p-i) nu(- p + q 2 +P ;p-^ D,{p- q ) 

-IPX 5 - p) KU - p) d,(p - 
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Neglecting would therefore be a suited approximation, especially since it is known to 
have very small impact on hadronic observables (although the vertex function is then, 
strictly speaking, no longer Poincare invariant). We will however follow another strategy 
here by symmetrizing our kernel in (q,p, P) <-» (— q, —p, P) for which we obtain 



-j[ri].,(^f±^-p)[ri] 



p + q + P 



;p-q) D w (p-q) 



p + q — P 



p + q — P 



p-qj [r 

q-p) [T J 



p + q + P 



p + q — P 



q-pj D^p-q) 
p-q) D n (p - q) 



p + q-P , . . 

nlru i ~ ; q ~ P ) D n (p - q) . 



(14) 



This is the only kernel we have found to analytically respect the axWTI for P 2 — > in the 
chiral limit for the general structure of the Bethe-Salpeter vertex given in Eq. (|12p . From 
a theoretical viewpoint it is less well grounded than the choice Eq. ( fl3l) : whereas ( fl3l) can 
be represented by proper Feynman diagrams the symmetrized version Eq. ( fT4l) does not 
share this property. It is however motivated from a pragmatic point of view since it allows 
one to check the influence of explicitly. We will therefore prefer the version Eq. ( fl3|) in 
this work. As a result we will see that neglecting contributions from H n introduces errors 
on the level of a view percent into observables such as the pion decay constant. For many 
applications it may therefore be possible to use the simpler kernel ([13]) . Both kernels have 
the correct charge conjugation properties and respect multiplicative renormalizability of 
the BSE by construction. 

The resulting diagrammatic expression of the Bethe-Salpeter equation is shown in Fig.0 
The sum over all contribution in Eq. f)14p is presented by the one-pion exchange diagram. 
Our approximation of the quark-DSE in Fig. by the one given in Fig. (jSJ) is justified 
provided the Bethe-Salpeter vertex functions of the pion with or without the pion inter- 
action term do not differ strongly. We will show in section IIIIBI that this is indeed the 
case. 

With the kernel (1141) our final expression for the quark-SDE is then given by 



s-\p) 



Z 2 S 1 (p) + g 2 C F Z 1F J -^-jL 7/i s (q) 7„ ^ 



k^kA Z(k 2 )T YM (k 2 ) 



k 2 



k 2 



-3 



d 4 q 



+r 7r 
+r 7r 



P + q \ or \r ( P + ( l 

— ;p-q)S{q)T 7T [ — ^— \q - p 

p + q 



2 

p + q 
p + q 



;q-p\ S(q)T w ( i^tl-p- ,, 



)p-q) S(q)T w ( p 



q-p) S(q) T n ( -^^-;p - q 



(15) 
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FIG. 6: The Bethe-Salpeter equation corresponding to the quark self energy of Fig. [5] 



with k = p — q. The corresponding expression for the BSE reads 

r*(p; P) = J k > p ) + *2£(p> p )) [S{k+)m P)S(k-)] sr (16) 

with the kernels Kjjf s and i^f^° r " given in Eqs. (jUJ) and (JHj). The Bethe-Salpeter vertex 
of the pion, i.e. the dressing functions E n , F n , G w and H v are normalized according to the 
condition 

2^P, = 3 Tr d j^ {rKfc, -P) ^ ( ^ P/2) ri(A:,P)g(A: - P/2) 

+K(k, -p)S(k + p/ 2 )r«*, p) ^/ /2) } 

+ /S/S^W( 9 ,-,)^^^ixiL(^) (17, 

with xi(^) P) — S(q + P/2)r{(/c, P)S(q — P/2). The condition is written for the momen- 
tum partitioning £ = 1/2 without loss of generality. The trace is over Dirac matrices and 
the conjugate vertex function T is defined as 

T(p, -P) = CT T (-p, -P)^- 1 (18) 

with the charge conjugation matrix C = —7274- The normalization condition ([IT]) guar- 
antees that the correspon ding residue of the pion pole in the four-point quark- ant iquark 
Green's function is unity |29j. The leptonic decay constant characterizing the pion cou- 
pling to the point axial field is then given by [30|] 

Hi, 

(19) 



f*=^T Tr * J -^f,{k-P)S{k + P/2) l5 PS{k-Pl2) 



P 2 =-M% 



where again the trace is over Dirac matrices. The numerical details needed for solving 
Eqs. ffl"5l) and ffT6l) simultaneously are discussed in subsection IIHAI We wish to emphasize 
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that the approximations made in the course of this section respect multiplicative renor- 
malizability. This is shown explicitly in appendix El Since also the axWTI is preserved 
we may hope that our approximation scheme captures essential properties of full QCD. 
In this context there are some interesting points to note: 

• A crudely simplified version of our interaction was considered by Gribov in his in- 
vestigations of quark confinement due to supercritical charges [3j. He conjectured 
that the inclusion of the pion back reaction to the quark induces important changes 
in the analytic structure of the quark propagator. We will come back to this point 
in subsection IIII Dl 

• We observe the quark self energy in Eq. (|T5|) to be infrared divergent in two di- 
mensions for any non-vanishing Bethe-Salpeter amplitude in the chiral limit. More 
precisely the axWTI requires r^(p, 0) = r l 7 5 5(p) / f n [B| and we have in good ap- 
proximation 

n . _ 3 [** 2 -'^+5^ (*) 



(27r) d V U J q 2 A 2 (q 2 ) + B(q 2 ) k 2 

for d dimensions. For any finite value of the dressing function B(p 2 ) at any momen- 
tum p, the integral of the pion contribution is infrared singular for d = 2 and would 
shift the value B(p 2 ) — > — oo, i.e. would give an infinite repulsive back reaction. 
Thus only the trivial solution B(p 2 ) = is left in two dimensions. This supports 
Coleman's theorem jiJ], which states that continuous symmetries in two dimensions 
cannot be broken dynamically. 

• Using the approximation ^(p, P) ~ r l/ ~f 5 B(p)/ f n \^ in the s-channel contribution 
of the pion to the quark- ant iquark scattering kernel (see Eq. ffTUl) ) a simpler choice 

for K Z™si<l,P\ p ) could be 

P) = -3[t j 75^A Ar(p - q) [r'Tb-y^U • (21) 

Jit Jit 

We consider this as an approximation to our construction for the kernel. 

• Finally we note that the inclusion of a back reaction by Goldstone bosons is possible 
for any case of dynamical symmetry breaking on an equal footing, e.g. it can similarly 
be applied to color-superconducting phases elaborated in [3l| for our approach. 

In the following subsection we will reconsider our approximation scheme from a different 
perspective. 

C. The l/iV c -expansion 

The approximation specified in the last subsection may also be viewed in the light 
of the l/iV c -expansion. The pion exchange contributions are then the (approximate) 
next-to-leading order correction to the leading one-gluon exchange of the pure Yang- 
Mills interaction. Such a view has been advocated already in model frameworks such 
as the NJL-model 0, [l(J In the following we discuss this aspect also for QCD in 



the framework of two-particle irreducible (2PI) effective actions, i.e. the CJT-formalism, 
which also allows for a derivation of the SDEs for propagators. The CJT action [32j for 



FIG. 7: Leading (sunset diagram on the left) and next-to-leading (right, the sum is over the number of gluon 
lines) order contribution to ^[S] in a model with effective one-gluon exchange. All signs and prefactors have 
been absorbed in the diagrams. 



truncations with effective one-gluon exchange is a functional of the quark propagator and 
given by 

T[S] = -TrLn (Z^SoS- 1 ) + Tr (l - Z 2 S 1 S) + T 2 [S] . (22) 

Here we already neglected the dependence on the expectation value of quark fields (ip) 
and (ip) as they are vanishing for the ground state. 

For the diagrammatic expansion of the 2PI functional r 2 [S'] we can then use a 1 /Re- 
ordering. The l/iV c -expansion of vacuum- vacuum contributions needed for l^fS] can 



nicely be arranged by the topology of the contributing diagrams [33J]. Since we do not 
consider pure gluonic contributions at order O(N^), the leading contribution is given by 
planar gluonic diagrams with a quark line as a boundary, being 0(N C ). The corresponding 
contribution to r 2 [S] is the left diagram of Fig. [7J To next-to-leading order, i.e. 0(1) , we 
have the topology of a cylinder with two quark lines as boundaries. Those contributions 
to ^[S] are subsumed in the right diagram of Fig. [7J 

Therefore the truncated quark-SDE to next-to- leading order in a l/iV c -expansion, given 
by = 0, turns out to be 

S-\p) = Z 2 S \p) + E<»(p) + E«(p) , (23) 
where E°(p) = J^( rambow ) (p) i s the rainbow self energy from Eqs. (EH) and 

*S>M = - /A«^fe f -? ;P -l) (24) 



(2*0 
with 

M tu , rs (q,p;P) = K™(q,p;P) 

+ J^K^(q f k ] P)S m (k + )M^ a (k fP] P)S b0 (k-) (25) 

is the next-to-leading order contribution. 

Since we do not want to perform the ladder resummation Eq. (1251) . we furthermore 
want to approximate M tu . rs (q,p; P) in order to reduce the numerical complexity. First 
we note that Eq. ( 1211) is still correct up to next-to- leading order in 1/N C when using 
the propagators from the rainbow-ladder approximation. The matrix M tu;rs (g,p; P) then 
contains Goldstone bosons in the s-channel. As we again presume the Goldstone bosons 
to be dominant, we will focus on the leading term in Eq. ([TO]) . Following our discussion 
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FIG. 8: The Bethe-Salpeter kernel K\ including next to leading order contributions in the pion interaction, cf. 
Eq. (|32[) . The dotted line in this context constitutes the resummed ladder being introduced in the l/./V c -expansion. 
In particular it includes the pseudoscalar channel. 



how to obtain the pion contribution to the quark self-energy given in Eq. (|15p from the 
s-channel contribution in Eq. ( flOl) . we can similarly motivate the same result for the self- 
energy as given in Eq. (T21]) : When using Eq. (ITU1) in Eq. (TS1]) , we are unable to find a suited 
kernel for the Bethe-Salpeter equation which respects chiral symmetry. The axWTI then 
suggests to use 

M tn . rs (q,p; P) 
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P ) p2 i Ml i^nlrs 

' 7T 


[-p; 


-P) + 




-?; 




[-p; 


P) 



(26) 

for the ladder resummed fully-amputated quark-antiquark scattering amplitude. This is 
again a pragmatic approach as we need to symmetrize the arguments in M tu . rs (q,p; P) 
in order to find a manageable Bethe-Salpeter kernel. We then recover the quark-SDE in 
Eq. (Tl5|) as the approximate next-to-leading order expression within the CJT-formalism. 

Strictly speaking we would then need to use the rainbow-ladder Bethe-Salpeter vertex 
functions in the expression for the self-energy, however to be consistent up to next-to- 
leading order in 1/N C we can also insert next-to leading order expressions. A simpler 
truncation would be to use the approximation given in Eq. (I2T]) . 

The Bethe-Salpeter kernel is then obtained by an additional functional derivative [34 

Ktu;rs = 7n ■ (27) 

0o rs 

Using our approximated next-to-leading order contribution of the self energy and neglect- 
ing the dependence of the Bethe-Salpeter amplitudes on the propagators then allows to 
justify the kernel given in Eq. ( tHj) . Essentially the kernel guarantees the pion to be 
massless in the chiral limit. 

Taking the l/A^-expansion more serious by performing the ladder summation results 
in a more complicated expression for the Bethe-Salpeter kernel. We need to evaluate the 
derivative in Eq. (T2"T|) . Also including the momentum in the multi-indices we abbreviate 
Eq. AM]) by 

vM - 



-S m Mi 



tu:rs 



We then get 



K 



(i) 

tu;rs 



M ts :ru + S x . 



xii, - M-tx;yu 
ODr.? 



(2c 



(29) 
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for the next-to-leading order contribution to the Bethe-Salpeter kernel. The leading order 
contribution is simply given by the effective one-gluon exchange. Evaluating the derivative 
by using Eq. fl25l) yields 

r Mt x -y U -^tx;ms-^rb;yu^bm Ki x - Tn M as -y U S na -\- 
0O rs 

Ktx;maS m aSbm~^~g ^ab;yu ■ (30) 

Now we can use Eq. ( 125]) to resum the right hand side of this equation. We then end up 
with 

r q -Mt X ;yu -^tx;ms-^rb;yuSbm -Mt X ;m-M aS ;y U S na . (^-Q 

0O rs 

Summarizing this analysis we obtain the next-to-leading order contribution to the Bethe- 
Salpeter kernel to be of the form 

■K-tu;rs M-ts;ur ■^tx;ms-^-rb]yuSbmS X y ~\~ M.i x - rn M. a s;yv,S naS xy j (*^) 

which corresponds to a resummed ladder in the t-channel and two contributions with two 
resummed ladders in the s-channel. The contributing diagrams are actually similar to 



those introduced in [121 ] . Having discussed these terms we can think of Eq. (1261) as an 
approximation of the first term of Eq. ( 1321) . which is identical to the one introduced in 
subsection III Bt A diagrammatic presentation of Eq. (1321) is shown in Fig. [HJ 



D. Gluonic interaction 



As discussed in previous subsections we need to specify two different components of the 
Yang-Mills part of the interaction: the gluon dressing function Z(p 2 ) and the vertex dress- 
ing Tym- The gluon dressing function has been calculated numerically from a truncated 



version of the Schwinger-Dyson equations for the ghost and gluon propagator [13l . 135 



The resulting dressing functions for the ghost and gluon propagator have been discussed 



and compared with corresponding lattice results in previous works, see e.g. |36j and refer- 
ences therein. In order to make this paper self-contained we use functional forms for the 
ghost and gluon dressing functions which represent the numerical solutions to sufficient 
accuracy 23]. These forms are given by 



2\ 



2 \ 2k / 

P \ I a mP 



G(p >) = ( *L_r72eW) (33) 



with the running coupling 



2 q(0) 4tt p 2 / 1 A 2 YM i 
"fltlP ) = \ , „9/a2 + IT Z3. , a2 ,„/„9/a2 x - 3 T2— > ( 34 ) 



\+p 2 Ia\ m Po P 2 + A 2 YM \]n(P 2 /A YM ) p 2 -A 2 



YM 



and the renormalization condition a(/z 2 ) = 0.968. The scale A Y m = 0.510 GeV is scheme 



dependent, here it corresponds to the momentum subtraction (MOM) scheme used in [13 
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FIG. 9: SDE solutions for the gluon dressing function Z(p 2 ) and the ghost dressing function G(p 2 ) adapted from 
ref. [35l ] and compared to the lattice results of [39l . |40| . 



351 ] . Via the analytic structure (13"3"j) of the gluon dressing function it can be related to 



the physical energy scale where positivity violations in the gluon propagator occur 
The infrared exponent k has been determined in an analytical infrared analysis (13 
k = (93 — vl20T)/98 ~ 0.595. The ultraviolet anomalous dimension of the ghost is 
given by 5 = —9N C /(4AN C — 8Nf) and related via 7 = —1 — 25 to the one of the gluon. 
Furthermore we have (3q = 4/7 m = (11N C — 2Nf)/3. Here we use N c = 3 and Nf = 2. 
Unquenching effects due to quark-loops in the gluon polarization are included in the 
expressions ( 1331 . 

We show together with lattice results for the unquenched (Nf = 2 + 1) gluon 



dressing function 3^] and the quenched 1 ghost dressing function [io| in Fig. For the 
gluon both results are in good agreement in the infrared momentum region (deviations 



are discussed in 36|). In the intermediate momentum region truncation artifacts in the 



gluon-SDE lead to an underestimation of the hump in the dressing function. In the ultra- 
violet the SDE-solutions show the correct logarithmic scaling as expected from resummed 
perturbation theory, whereas the lattice data suffer from artifacts due to the finite lat- 
tice spacing. The ghost dressing function show slight deviations in the infrared, which 



are discussed in detail in [36|. In general we believe that both solutions are in sufficient 
agreement such that the expressions (13"3l serve as a reliable input for the quark SDE. 

For the second part of the Yang-Mills interaction, the vertex dressing T Y m, we use a 
procedure first suggested in [25] and further explored in 26[: we devise a functional form 
for the vertex dressing with parameters fixed by the requirement that lattice data for 
the quark propagator are reproduced by solutions of the quark-SDE using the complete 
interaction. The functional form we employ for this procedure is given by 



7i 



,r yM (^,/i 2 ) 



7,rMA; 2 )r : 



(35) 



1 Lattice calculations of the unquenched ghost dressing function are not yet available. 
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with the components 



^- Mwm^y (36) 



k 2 ~ 
,(/A,r) = J k 2 + A 2 YR G(k 2 ^ 2 ) G(CV) Z^ 2 ) h [\n(k 2 /A 2 YM + r)] 1+5 (37) 

o o o a(-W) + k 2 /A 2 nrn 

T 3 (k 2 , ^ 2 ) = Z 2 ^ 2 ) \[ h2 ' 2 , (38) 

1 + K I^QCD 

where r = e — 1 acts as a convenient infrared cutoff for the logarithms. As mentioned 
in section III AI other tensor structures in the Yang-Mills part of the quark-gluon vertex 
might be relevant, an investigation of those is however beyond the scope of this work. We 
will come back to this point at the end of section III Dl In the following we discuss the 
choice of our ansatz ( |35j) step by step and compare it to the one used in the quenched 



calculation of ref. 26 



It is well known that the effective interaction in the quark-SDE has to approach the 



running coupling in the ultraviolet momentum regime [41 



i.e. 



f -V^(A: 2 )r 1 (P)r 2 (A; 2 )r3(P) - J 7 " (39) 
4tt z 2 Zz \n(k 2 /A 2 QCD ) 

Up to constants this UV-behavior is represented by Ti with the scheme- dependent scale 
^■qcd- The anomalous dimension j m of the quark propagator is given by 7 m = 12/(33 — 
2iVy), where we use Nf = 2 in our calculations. Furthermore note that Z 2 = Z 3 Zip in 
Landau gauge, where Z 3 is the renormalization factor of the ghost fields. 

The product Z(k 2 )T 2 (k 2 )T 3 (k 2 ) goes to a constant for large momenta; the scale A Y m 
is determined from the SDE-solutions for the ghost and gluon propagator fl33|) . The 
coefficient h is fixed such that the limit (139]) is satisfied. The renormalization group 
invariant G(( 2 , fi 2 )Z 3 (fi 2 ) with the arbitrary but fixed scale ( is introduced to impose the 
correct cutoff- and renormalization point dependencies of the effective interaction in the 
quark-SDE. In our calculations we use ( = 2.9 GeV, other choices do not affect the result. 

The infrared behavior of our interaction is dominated by T 2 , which is diverging for 
p 2 — > whereas Ti and I^ go to constants. The degree of divergence is given by 
G{k 2 ) <l ' k 2+^2 ~ (/c 2 )~ K ~ 1 / 2 . This behavior has been derived in ref. |24j from an an- 
alytic, selfconsistent analysis of the full SDE for the quark-gluon vertex given in Fig. [lj 
We wish to emphasize, however, that the precise infrared behavior of Tym is not impor- 
tant for the quark-SDE: this equation is dominated by intermediate loop momenta. One 



could equally well work with lower degrees of divergence [26J or even with an infrared 
finite vertex (25| . 

Finally, from the Slavnov- Taylor identity of the quark-gluon vertex one can see that the 
vertex is also proportional to the quark wave function Zf. This dependence is taken care 
of by T3, which mimics the momentum dependence of Zf. The extra factor Z 2 is vital for 
ensuring multiplicative renormalizability of the quark-SDE. The dependence of this part 
of the vertex on the quark mass is expressed in terms of the function 

° (M) = 1 + a 2 M(( 2 )/A QCD + a,M 2 (C)/^ Q cD ^ 
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FIG. 10: SDE results compared to lattice data of the unquenched quark mass function M(p) and the wave function 

z f (p) m. 



h 


A rM (GeV) 


A QCD (GeV) 


ai 


a 2 


a3 


0.99(1) 


0.51(5) 


0.65(5) 


5.22 (1) 


5.06 (1) 


-9.06 (1) 



TABLE I: Parameters used in the vertex model, Eqs. (|35l40p . 



In order to preserve multiplicative renormalizability of the quark-SDE, it is important 
that the scale ( is not correlated with the renormalization point. Instead it should be a 
fixed scale. In our calculations we use the same scale ( = 2.9 GeV as also present in T 2 . 

To fit the parameters of the Yang-Mills part of our interaction we use the following 
procedure: We first determine the values of h,Aq C D and a(M) at a fixed current quark 
mass such that (i) the correct ultraviolet behavior ( 1391) of the running coupling is repro- 
duced and (ii) the unquenched lattice quark propagators from ref. [15J are reproduced. 
The Yang-Mills scale Aym is taken from the SDE- results parameterized in Eqs. This 
procedure is repeated for the four different current quark masses of ref. [HI , which corre- 
spond to Af (2.9 GeV) = 44, 65, 85, 106 MeV. We then fit Eq. gO} to the results for a(M) 
and determine a\, a-i and a^. The results for the parameters are given in Tab. [H 

The numerical results of our fitting procedure are shown in Fig. [TUJ For the quark mass 
function M(p 2 ) we find excellent agreement with the lattice data in both, the infrared 
and ultraviolet momentum regions. For the wave function Zf(p 2 ) the agreement is only 
slightly less convincing. Particularly interesting is the small spread of the wave functions 
for different quark masses, which is reproduced by the SDE results. This is in contrast 
to earlier investigations of the quenched quark propagator, where the spread in the SDE 
results has been too large 25, 26[. We consider the improved result here as an indication 
that our quark-gluon interaction is more realistic than in previously used models. In 
particular we attribute this to the fact that the hadronic contributions in the quark-gluon 
interaction generate all twelve possible Dirac structures of this interaction (only eight of 



26 



these are independent in Landau gauge), whereas the quenched investigations [2 
mainly worked with the 7^-structure only. This is a clear hint that one should include 
further structures in the modeling of the Yang-Mills part Tym of the interaction and is left 
for future work. Our solutions for the quark mass function M(p 2 ) and the wave function 
Zj{p 2 ) for physical up/down quarks are discussed in section [Til Bl 
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FIG. 11: The momentum routing in the SDE and BSE. 

III. NUMERICAL PROCEDURE AND RESULTS FOR QUARK AND PION 

A. Numerical Procedure 



Before presenting results we need to detail our numerical procedures. Our choice of 
the momentum routing in the SDE and BSE is specified in Fig. [Hj In the right diagram 
P denotes the total momentum of the meson and q + = q + P/2, g_ = q — P/2 and 
p + = p + P/2, p- = p — P/2 are the momenta of the internal and external quarks. Since 
in the rest frame of the meson we have P M = (0, 0, 0, iM), where M is the mass of the 
bound state, we need to know the quark propagators in the BSE for complex momenta. 
In addition, for the pion exchange diagram we need to know the pion Bethe-Salpeter 
amplitude for complex relative momenta. Whereas solving the quark- SDE in the complex 
plane is feasible and a standard procedure by now (see e.g. [IH, |H, 43]), the solution of 
the BSE for complex momenta has not yet been performed to our knowledge. In view 
of the considerable numerical complexity of the coupled system of SDE and BSE we will 
not attempt such a solution within the context of this work. Instead we resort to the 
well explored 'absolute value approximation', which amounts to replacing all momentum 
arguments in internal quark propagators in the BSE by their absolute value. In addition, 
we replace the arguments of the pion amplitudes in the meson exchange diagrams by 
their real parts. This approximation has the merit of making the kernel independent of 
P^, which in turn greatly simplifies the Bethe-Salpeter norm (TTTj) . We will attempt a 
complete solution of the problem in the complex plane in future work. 



To estimate the quality of this approximation we employed the (quenched) model of [43 
and solved it once for complex momenta and once in the absolute value approximation 
for a physical pion. Using similar current quark masses the resulting difference in the 
pion mass is less than 1 MeV. The error in /„. however is larger: we obtain /„. = 94 
MeV with complex momenta but only f n = 82 MeV in the absolute value approximation. 
Roughly 10 MeV of this difference is due to the absolute value approximation in the norm 
integral (|17|) : the remaining 2 MeV are accounted for by shape distortions in the F, G 
and H amplitudes of the pion Bethe-Salpeter amplitude ffT2l . We expect similar errors 
when using the quenched version of our interaction, i.e. only taking into account the 
gluon exchange part of the interaction given by Eq. (135]) . In the unquenched calculations 
the error in the norm integral is backfeeding into the coupled system by the overall 
normalization of the pion contribution of the interaction. We therefore expect somewhat 
larger errors in both, the pion mass and the resulting pion decay constant. Nevertheless 
we consider our results to be accurate and meaningful in their qualitative features. 

Further details concerning the numerical techniques needed to solve the quark-SDE are 
given in 13, lH; corresponding techniques for the BSE are described in jol 141 . 
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FIG. 12: In the upper panel we show the quenched and unquenched (Nf — 2) quark mass function (left diagrams) 
and wavefunction (right diagrams) for physical up/down quarks with M(2.9 GeV) = lOMeV. In the lower panel 
we compare the SDE/BSE-solutions for a heavier quark of M(2.9 GeV) = 44 MeV to the quenched and unquenched 
lattice results of Bowman et al. [l5l |. 



B. Quark propagator and pion properties 

Our results for the physical up/down quark propagator are shown in Fig. [12j We 
compare the unquenched (Nf = 2) solutions from the SDE/BSE system with partially 
quenched solutions. Here, 'partial quenching' means that we still include the quark loop 
effects in the gluon propagator, implicitly present in Eqs. (1331) . but switch off the hadronic 
unquenching effects from the pion part of the interaction. In both calculations the renor- 
malized current quark mass is chosen such that the pion mass is M n = 136 MeV. In an 
MS-scheme with Ajfg = 225 MeV this corresponds to renormalized up/down current 
quark masses of 

m^ =2 (2 GeV) = 3.5 MeV rn"f=°(2 GeV) = 2.8 MeV. (41) 

Note that the unquenched current quark mass is larger than the quenched one. In [13] 
a small effect in the opposite direction has been found from quark loops in the gluon 
polarization alone. The stronger increase found here reverses this tendency and indicates 
that unquenching effects in the quark-gluon vertex are stronger than those in the gluon 
polarization. 

The quark mass function in the upper left diagram of Fig. [12] also shows a sizable 
reduction of the dynamically generated mass M(p 2 ) in the infrared when the theory is 
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unquenched. These effect is greatly diminished at p = 2 GeV and vanishes for momenta 
larger than p = 3 GeV. The corresponding effect in the quark wave function Zf(p 2 ) 
is much smaller and only visible for momenta lower than p = 0.5 GeV. This finding 
agrees qualitatively with results from quenched and unquenched (Nf = 2 + 1) lattice 
calculations shown in the lower panel of Fig. [T5J Clearly, both approaches underpin the 
physical intuition: in general, quark loop effects are screening and compensate part of the 
antiscreening gluonic interaction. This leads to a reduction of the interaction strength in 
the gluon propagator 13, [39|. Here it also reduces the strength of the quark-gluon vertex 
and consequently entails (further) reduced quark masses. The effects for the quark mass 
function seen in the SDE/BSE approach are larger than in the lattice calculations. This 
may or may not be an artifact of the absolute value approximation used in the numerics. 
In any case, it seems save to conclude that unquenching effects from the pion back reaction 
(i.e. in the quark-gluon vertex) are considerably stronger than the quark- loop effects in 
the gluon polarization investigated in refs. [H, 14 



The corresponding Bethe-Salpeter amplitudes of the pion are shown in Fig. [T3J Sim- 
ilar to the previous findings in the quenched approximation 42| we find a hierarchy in 
magnitude of these amplitudes: clearly E n is most important in size followed by F n , 
and H n . This hierarchy is not changed by unquenching, although we note a considerable 
increase in size of G n and H v for Nf = 2 compared to the quenched case. Despite these 
differences both sets of amplitudes are similar enough to justify the approximation made 
in subsection III Bl In the infrared, E n , F n , and go to constants. In the ultravi- 
olet E n and F n fall like 1/p 2 times a logarithm, whereas G n and H n are proportional to 
1/p 4 times a logarithm. This behavior is not changed when the system is unquenched. 
The relative importance of the Bethe-Salpeter amplitudes can be inferred from table [TT1 
We show results for the quark mass function at zero momentum, the pion mass and the 
contributions to the pion decay constant when part of the amplitudes are omitted. On 
a quantitative level it is certainly important to include all four amplitudes in the calcu- 
lation. If, however, one is satisfied with qualitative results, even a calculation involving 
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TABLE II: The influence of the four components of the Bethe-Salpeter amplitude of the pion on masses and decay 
constants. All results are for the unquenched interaction; the units are in MeV. 



only the ^-amplitude would give satisfactory results. 

Using all four Bethe-Salpeter amplitudes, the resulting values for the pion decay con- 
stant calculated from (fl9l) in the absolute value approximation and for a physical value 
of the current quark masses are 

f"f= 2 = 49 MeV f? f= ° = 82 MeV. (42) 

Clearly, unquenching reduces the value of f n . However, the magnitude of this reduction 
is too large. Also both values are too low. One would expect the unquenched value close 
to the physical point f^ xp = 93 MeV and a somewhat larger number in the quenched 
case. A reduction of the unquenched pion decay constant was also found in recent lat- 
tice calculations with quenched [Hf and unquenched |45| twisted mass fermions. They 
obtained 

fnSL = 86(1) MeV f^ZL = 102(3) MeV. (43) 

for our definition of the pion decay constant and in the chiral limit. In the unquenched 
simulation the pion decay constant at the physical point was used to determine the phys- 
ical scale. The uncertainty of this scale, at least in quenched calculations, is usually 
expected to be on the order of 10%. 

As discussed in section IIII Al at least part of the deviations of ( H2|) from (143]) can be 
attributed to the absolute value approximation. However, it is by no means clear, whether 
a calculation in the complex plane will bring /„- all the way up to the physical value. An 
indication that this may well not be the case is a comparison with /„. from a modified 



Pagels-Stokar approximation [46|, |47 



a _ N c „ [ , 2 2 M(q 2 )Z f (q 2 ) / „ 2 . q z dM{q 2 '' 



a,* # * jffwm \ M(q) ~ ■ (44) 

This approximation is not affected by the absolute value problem, but in turn incorporates 
only the effects of the leading pion Bethe-Salpeter amplitude in the chiral limit. Thus 
the systematic errors are different. By comparison of (jUj) with the full result in model cal- 



culations [42j, |43j one finds that the approximation (|44l) should lead to an underestimation 



of the full f n of the order of 10%. In our case we obtain 

/5>= 2 = 65 MeV tf/ P =° = 82 MeV. (45) 
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Compared to (j42j) we find a considerable increase in the unquenched value and no change 
for the quenched case. Within our quenched interaction we can also perform a full cal- 
culation in the complex plane using the techniques described in We then find the 
value 

= 89 MeV, (46) 

which is indeed roughly 10% larger than the result ( l42|) of the absolute value approxi- 
mation and also 10% larger than the result ( l4"5l) from the Pagels-Stokar approximation. 

From this one would expect to find a value of roughly f Wi compiex ~ 70 — 75 MeV from an 
unquenched calculation in the complex plane. Clearly both, the calculated result ( 1461) for 
the quenched case and the conjectured result for the unquenched case fall short of fj4~3l) 
by almost twenty percent. 

A possible source of these shortcomings are errors in the overall scale of our calculation. 
This scale was set by comparison with the lattice results for the gluon 3^] and quark 



propagator [15j. Indeed, already in [25j and again in [26| it was found that the pion 
decay constant calculated from an interaction fitted to quenched lattice data yields an 
underestimation of by about thirty percent. This agrees with our findings above and 
supports the argument that there may be problems with the lattice scale. We also note 
that a rescaling of the lattice units would lead to larger constituent quark masses and 
in addition reduce the effect of the Goldstone contribution in the selfenergy. This would 
improve the agreement of the results shown in Fig. (fl~2]) . 

C. Gell— Mann-Oakes-Renner relation and chiral condensate 

The Gell-Mann-Oakes-Renner relation 

Ml = [mM + m d Gu)]<^>V) (47) 

can be derived from the inhomogeneous Bethe-Salpeter equations for the axialvector and 
pseudoscalar currents and the axialvector Ward-Takahashi identity [7| . Since the left-hand 
side of this equation only contains physical observables, the dependence of the renormal- 
ization point of the current quark masses m u (fi) and md(/i) and the chiral condensate 
have to cancel each other. The individual dependence of the pion mass and the 
pion decay constant on the current quark mass is shown in the upper panel of Fig. O 
Clearly, one observes the characteristic square-root behavior of the pion mass that be- 
comes massless in the limit m u /^ — > 0. No qualitative changes can be seen by naked eye 
between the quenched and unquenched calculation. The situation is somewhat different 
for the pion decay constant. Here an almost linear behavior for sizable quark masses can 
be seen, with small deviations towards the chiral limit. Whether the different curvature 
for m M /rf — > can be explained by the presence of chiral logs remains to be explored in a 
more refined calculation. In the lower panel of Fig. [TD we plot M% as a function of the 
current quark mass. Clearly, the dependence is linear for a wide range of current quark 
masses. The deviations from the linear fit for small quark masses may well be an artifact 
of the increasing numerical uncertainties close to the chiral limit. The deviations for larger 
quark masses show the onset of quantitative corrections to the Gell-Mann-Oakes-Renner 
relation. It is interesting to note, however, that these corrections play no role up to pion 
masses of M n = 400 MeV. A similar result has been found in a recent lattice study using 



Wilson quarks and fine lattices 48 
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FIG. 14: The quenched and unquenched (TV/ = 2) pion mass and decay constant as a function of the renormalized 
quark mass m(fj,) at \i = 2 GeV (MOM-scheme). The straight lines in the lower panel show the validity of the 
Gell-Mann-Oakes-Renner relation. 



The chiral condensate can be extracted directly from the quark propagator via 

(^iW), (48) 

where the trace is over Dirac indices, Sq is the quark propagator in the chiral limit and Z m 
is the quark mass renormalization factor that can also be determined from the quark-SDE. 
We obtain 

(«^ =2 (2 GeV) = (-240 MeV) 3 (**)^ /=0 (2 GeV) = (-283 MeV) 3 , (49) 

Ivl ij o 

which shows a sizable unquenching effect in the condensate. Whether this effect stays as 
large when the absolute value approximation is given up remains to be seen. Nevertheless 
it is instructive to use these values as input into the Gell-Mann-Oakes-Renner relation. 
This allows us to extract the pion decay constant once more from an independent source. 
We obtain 

fnfMOR = 74(3) MeV f%£ OR = 83(3) MeV, (50) 
for the absolute value approximation and 

fH,GMOR,com P le X = 85(3) MeV, (51) 
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FIG. 15: The absolute value of the Schwinger function as(t) of the quark propagator. 



for the quenched calculation using complex momenta. The error subsumes numerical 
errors as well as errors in fitting M 2 (m), which has been performed by a least squares 
fit in the region from m(2GeV) = 3 — 40MeV. By comparison with (|46|) we note that 
the GMOR- relation is satisfied at the 4 percent level for the central value in (l5"Tj) . The 
difference of the corresponding value for the unquenched case and the result Eq. (T4"2"j) shows 
once more that the absolute value approximation leads to a large error in the unquenched 



Finally we take a look at the analytic structure of the quark propagator in the complex 
momentum plane. We wish to stress that the results of this section have to be treated 
with caution, since the absolute value approximation is likely not to be reliable in this 
respect. Nevertheless, it may be instructive to have a first glance in the complex plane. 
To this end we determine the Schwinger-f unction 



where crsip 2 ) = B{p 2 ) / {p 2 A{p 2 ) 2 + B{p 2 ) 2 ) is the scalar part of the dressed quark propa- 
gator. (This method has a long history, see e.g. [23J and references therein). According 
to the Osterwalder-Schrader axioms of Euclidean field theory this function has to 
be positive to allow for asymptotic quark states in the physical sector of the state space 
of QCD. Conversely, positivity violations in the Schwinger function show that the corre- 
sponding asymptotic states (if present) belong to the unphysical part of the state space. In 
Fig. [15] we show results for the Schwinger function of the quenched and unquenched quark 
propagator for the physical up/down quark of Fig. [T21 Both functions are qualitatively 
similar. An excellent fit to the Schwinger function is obtained using the form [50(] 



D. Analytical properties of the quark propagator 




(52) 



*s{t)\ 



6 exp(-&it) cos(b 2 t + b 3 ) 



(53) 
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which corresponds to a pair of complex conjugate poles of the propagator in the time-like 
momentum plane. These poles correspond to a 'quark mass' given by m = b\ ± ibi- In 
our quenched case this amounts to m = 505(10) ±i 360(10) MeV; the unquenched quark 
mass is m = 351(10) ± i 254(10) MeV. The mismatch of the fit with a(t) at small times 
shows the presence of additional analytic structure either in form of a cut along the real 
time-like momentum axis or in form of additional singularities further away from the 
origin. If this behavior is stable also when the absolute value approximation is overcome 
it contradicts the Gribov scenario given in j3[. However, as emphasized already above, 
definite conclusions can only be drawn from a complete calculation including the full 
complex momentum structure involved in the coupled system of the quark-SDE and the 
pion-BSE. 



IV. SUMMARY AND OUTLOOK 



In this work we investigated the pion back reaction on the quark propagator and result- 
ing pion properties. To the best of our knowledge this is the first calculation of this kind in 
a non-perturbative continuum approach to QCD. We isolated the pion contribution to the 
quark-gluon vertex and identified an approximation to the quark-SDE which allows for 
the construction of a Bethe-Salpeter kernel in agreement with the axial Ward-Takahashi 
identity. This setup is powerful enough to analytically verify Coleman's theorem in 1+1 
dimensions: chiral symmetry cannot be broken spontaneously there. Our main interest, 
however, is QCD3 + i, where we evaluated the back reaction effects numerically. Here we 
had to resort to an absolute value approximation, which is reliable as concerns the pion 
mass but problematic for the pion decay constant. 

As a result we have obtained considerable unquenching effects in the quark mass function 
at small momenta. The screening effect of the pion interaction reduces the quark mass 
M(0) from M N f =0 (0) = 477 MeV to M N f =0 (0) = 371 MeV. This effect becomes smaller 
for larger momenta and vanishes above p = 3 GeV in correspondence with the width of 
the Bethe-Salpeter amplitude of the pion. Almost no corresponding effect in the quark 
wave function is seen. 

The screening effect of the pion interaction also reduces the value of the chiral con- 
densate. The Gell-Mann-Oakes-Renner relation is satisfied with an accuracy less than 4 
%. We obtain an almost linear relation between the squared pion mass and the current 
quark mass for pion masses up to 400 MeV with only slight deviations up to 500 MeV. An 
unsolved problem is posed by our low value of the pion decay constant. Certainly, part of 
the problem is the absolute value approximation mainly used in this work. However, we 
also argued that the overall scale of our calculation, which is based on the lattice results 
of [11,0, may have to be corrected by about 20 %. 



Our findings are certainly reliable on a qualitative basis. To also obtain reliable quanti- 
tative results we need to go beyond the absolute value approximation for the arguments 
in the internal quark propagators and pion amplitudes. Work in this direction is well 
under way. We then hope to better understand the problem of the very low pion de- 
cay constant obtained in this work. Also, working in the complex momentum plane is 
mandatory to investigate possible unquenching effects in the analytic structure of the 
quark propagator conjectured by Gribov. Finally, a complex treatment of the system is 
mandatory to explore decay thresholds in vector, axialvector and the scalar meson sec- 
tors. An investigation of this mesons should be simplified by use of the Maximum Entropy 



25 



Method 51 . 
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APPENDIX A: MULTIPLICATIVE RENORMALIZ ABILITY OF SDE AND BSE 

Here we show explicitly that our approximation scheme for the Schwinger-Dyson and 
Bethe-Salpeter equations preserves multiplicative renormalizability (MR). To this end we 
rewrite (115j) and (fl6l) in symbolic notation and make explicit the dependence of the various 
quantities on the renormalization point fi. We begin with the SDE: 

S-\p,fi 2 ) = Z 2 ( f i 2 )S \p)+g 2 (fi 2 )Z 1F ( f i 2 ) J S(q,ii 2 )Z(k 2 ^ 2 )T YM (k 2 ^ 2 )L YM ( V , qi k) 

-3 J [T*(p, q, fi 2 ) S(q, /i 2 ) Hp, q, fi 2 ) L n (p, q, k)] , (Al) 

where the functions L YM and L w subsume all \i- independent quantities in this equation. 
Since 

S-\p, /x 2 ) = ip-yAtf, /x 2 ) + B(p 2 , /i 2 ) (A2) 
and we have the relations 

A{p 2 ^ 2 ) = Z 2 ^ 2 ,A 2 )A (p 2 ,A 2 ) 

B{p 2 ^ 2 ) = Z 2 ^ 2 ,A 2 )B (p 2 ,A 2 ) (A3) 

between the renormalized quantities A, B and the unrenormalized, cutoff (A) dependent 
quantities Aq, Bq, we note that the //-dependence of the left hand side of Eq. (1A1I) is given 
by Z 2 {ij j 2 )- This is also trivially true for the tree- level term Z 2 (/j, 2 ) Sq 1 (p). It follows then 
that the whole equation is multiplicatively renormalizable if and only if each term of the 
self energy is also proportional to Z 2 (fi 2 ). To show this we need the relations 

Z(k 2 ,fi 2 )Zj^ 2 ,A 2 ) = Z (k 2 ,A 2 ) 

G(k 2 ,fi 2 )Z 3 (^,A 2 ) = G (k 2 ,A 2 ) 

g(f, 2 )Z g (f, 2 ,A 2 ) = g (A 2 ) (A4) 

between the renormalized gluon and ghost dressing functions Z, G and their unrenor- 
malized counterparts Z , Gq and a similar relation for the coupling g. With the help of 
(IA4j) we extract the renormalization point dependence of the Yang-Mills interaction from 
Eqs.((35D-(l3HD: 

T YM (k 2 ^ 2 )-^- (A5) 

The renormalization point dependence of the Yang-Mills part of the self energy is then 
given by 

Z YM (p 2 , /i 2 ) . (A6) 

ZW)Z 3 ^)Z 3 (^) 



2G 

Taking into account the Slavnov- Taylor identities 

Zip = ZgZ^zJ , Z\ = Z g Z^zJ (A7) 

we arrive at 

S™(pV)~§^, (A8) 

which contains the ghost-gluon vertex renormalization constant Z\. In Landau gauge 
the ghost-gluon vertex is finite and consequently this quantity can always be chosen to 
equal one by a suitable renormalization condition 2 . Thus we arrive at the desired result 

The other contribution to the quark self energy trivially gives the same result provided 

T*(p, q ^ 2 )~Z 2 (^). (A9) 

This dependence agrees with the one required in the norm integral, Eq. ( JT7j) . and the 
expression for /„., Eq. ffl9|) and therefore completes the proof of MR for the SDE. 
In the same symbolic notation the Bethe-Salpeter equation for the pion can be written 

as 



r w (p,P) = <?V) J r 7r (g, P, /i 2 ) S(q-, fj 2 ) S(q+, /i 2 ) Z\p Z (k 2 , /i 2 ) Yyuik 2 , /i 2 ) Jym(p, Q, P) 
+ J 1^(9, P, // 2 ) 5( 9 _,^ a ) %+,A/ 2 ) [1^(9, p,P,// 2 )] 2 J ff (p, g ,P) (A10) 

where the functions Jym and again subsume all //-independent quantities. With the 
help of the relations given above it is easy to verify that all terms on the left and right hand 
side of this equation are proportional to which leaves the equation multiplicatively 
renormalizable. 
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